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We use bosonization to derive the effective field theory that properly describes ferromagnetic transi- 
tion in one-dimensional itinerant electron systems. The resultant theory is shown to have dynamical 
exponent z = 2 at tree level and upper critical dimension d c — 2. Thus one dimension is below 
the upper critical dimension of the theory, and the critical behavior of the transition is controlled 
by an interacting fixed point, which we study via e expansion. Comparisons will be made with the 
Hertz-Millis theory, which describes the ferromagnetic transition in higher dimensions. 
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Ferromagnetic transitions in itinerant electron systems 
are among the very first examples of quantum phase 
transitions studied theoretically [1]. In the approach 
pioneered by Hertz, one decouples the electron-electron 
interaction using Hubbard-Stratonovish transformation, 
integrates out the fermionic degrees of freedom, and ar- 
rives at a Ginsburg-Landau- Wilson like free energy func- 
tional that involves the ferromagnetic order parameter 
only, which are bosonic degrees of freedom. The quan- 
tum nature of the theory lies in the (imaginary) time de- 
pendence of the order parameter. This effective bosonic 
theory, known as the Hertz-Millis theory, was argued to 
have upper critical dimension d c = 1, thus the critical 
behavior of the transition is expected to be mean-field 
like for both d = 2 and d = 3 [1,2]. It has been pointed 
out recently, however, that the procedure of integrating 
out gapless fermions may lead to subtle singularities in 
the expansion of the resultant bosonic free energy func- 
tional in terms of the order parameter, or its gradients 
[3]; such singularities may invalidate the power-counting 
analysis of Hertz and Millis, and change the critical be- 
havior of the transition [3]. The nature of the transition 
in two- and three-dimensional itinerant electron systems 
is currently under extensive theoretical study [3] . 

Comparatively speaking, much less attention has been 
devoted to the possibility of ferromagnetic transition in 
one-dimensional (ID) electron systems, until recently. 
This is in part due to a theorem of Lieb and Mattis [4] , 
which states that the ground state is a singlet for certain 
classes of one-dimensional models with spin-independent 
interactions, thus rules out the possibility of ferromag- 
netism in these models. The existence of ferromagnetic 
ground states in ID models (that are not dictated by 
the Lieb-Mattis theorem) was only established recently 
through numerical work [5], in which it was also found 
that the ferromagnetic transition is second-order, thus 
there is a quantum critical point. On the experimental 
side, the so-called "0.7(2e 2 )//i" anomaly (or 0.7 anomaly) 
in the density-dependence of the conductance of a ID 
electron wire has attracted considerable attention; one 
of the possible interpretations of this anomaly is sponta- 



neous magnetization of the electrons, which would give 
rise to plateau-like behavior near 0.5(2e 2 )//i [6]. This 
interpretation receives further support from the observa- 
tion of similar behavior in the presence of magnetic field 
near crossings of subbands with opposite spin polariza- 
tion [7]; in this case the corresponding interpretation is 
spontaneous pseudospin magnetization of the electrons, 
where the pseudospin index is actually the subband in- 
dex. Inspired by these developments, a description of the 
ferromagnetic phase of a ID electron liquid is developed 
in Ref . [8] , although the critical behavior of the ferromag- 
netic transition was not studied in that work. 

In this paper we study the critical behavior of ferro- 
magnetic transition in ID. We derive the effective bosonic 
field theory that describes the transition using Abelian 
bosonization, and use renormalization group to study the 
critical behavior. The bosonization method, which is spe- 
cific to ID systems and extremely powerful, allows one 
to derive the effective bosonic theory in terms of the fer- 
romagnetic order parameter without having to integrate 
out gapless fermions; thus the resultant theory does not 
suffer from singularities that the Hertz-Millis theory may 
encounter. A straightforward power-counting analysis of 
the theory indicates that the dynamical exponent z — 2 
at the transition (at tree- level), thus the upper critical di- 
mension is d c = 4 — z — 2. As a result the theory is below 
its upper critical dimension in ID, thus the critical behav- 
ior of the transition is controlled by an interacting fixed 
point. We study the critical behavior using momentum 
shell renormalization group combined with e expansion 
near d c = 2, at zero as well as low temperature. 

Model and Bosonization — Consider the follow- 
ing Hubbard-like Hamiltonian describing interacting elec- 
trons in ID: 

H = % c L c j> + U Y1 n iT n H + Vij" m*n j<Tf , 

ijcr i ijcro-' 

(i) 

where tij are the single electron hopping matrix elements, 
U > is an onsite repulsion, and V™ represents (possi- 
bly spin-dependent) further neighbor repulsion; depend- 
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ing on the spin dependence of V the system may possess 
either full Heisenberg (or 0(3)) symmetry, or just Ising 
(or Z(2)xU(l), the latter is responsible for conservation 
of S z ot ) symmetry. It has been shown [5] that in the pres- 
ence of further neighbor hopping (ty with \i — j\ > 1), 
the Lieb-Mattis theorem no longer applies and a ferro- 
magnetic ground state is stabilized for large enough U , 
when V = 0. One can also stabilize the ferromagnetic 
phase by having V terms that are spin-dependent and 
ferromagnetic; this possibility has been considered in the 
context of atoms trapped in ID optical lattices [9]. 

One of the most powerful methods of tackling such ID 
models is Abelian bosonization [10]. In this scheme, one 
takes advantage of the fact that in ID, all particle-hole 
excitations can be generated by electron density and cur- 
rent operators which satisfy bosonic commutation rela- 
tions, and expresses both kinetic energy and interaction 
terms of H in terms of the electron density and current 
operators. If one keeps terms that are quadratic and with 
least number of gradients in the density and current op- 
erators, one arrives at the familiar Luttinger liquid (LL) 
Hamiltonian, which describe decoupled spin and charge 
excitations of the paramagnetic phase [10]: 

Hll = H c + H s ; (2) 
Hc = h I dx ^ v JcK{x) + VNc(d x Mx)) 2 }; (3) 
Hs = hj dx ^ v J^ x ) + VNs(dAs(x)?]- (4) 

Here 4> c (x) and 4>s{x) are the charge and spin fields re- 
lated to the charge and spin densities of the system: 

p(x) = -d x (t> c {x), S z (x) = ^-d x (j) s {x); (5) 

while II Q are their conjugate fields satisfying 

[4> a {x),Jl al {x')] = i6 aa >5(x - x'), (6) 

with a being c or s. Physically U a (x) represents local 
charge or spin current. Clearly the velocity parameters 
vmc and vn s parametrize the energy cost of charge and 
spin density fluctuations respectively, and are thus pro- 
portional to the inverse charge and spin susceptibilities, 
while v,jc and v,j s are proportional to the charge and spin 
stiffness of the system respectively, as they measure the 
energy cost of charge and spin current fluctuations. 

We emphasize that Hll is an approximation of the 
original electron Hamiltonian Eq. (1) with generic sin- 
gle electron dispersion relation and two-body interaction. 
For example, a nonlinear term of the form J dx cos(\/80s) 
that describes back scattering of electrons with opposite 
spins is neglected here. As pointed by Haldane [11], non- 
linearity in single-electron dispersion gives rise to terms 
beyond quadratic order in n and <f>, which represent in- 
teractions among the bosons. Also the non-locality of 



electron-electron interaction (due to V for example) leads 
to non-trivial wave- vector dependence in Fourier space, 
which gives rise to terms that are quadratic in <j> but 
involve higher gradients. These terms, however, are ir- 
relevant in the renormalization group sense, at the Lut- 
tinger liquid fixed point described by Hll (Eq. (2)); they 
scale to zero in the long- wave length and low-energy limit; 
their physical effect is to rcnormalize the parameters of 
Hll- Thus the long- wave length, low-energy properties 
of the system arc well described by Hll, albeit with 
renormalized parameters. This is the essence of the Lut- 
tinger liquid theory [11,10]. What we are going to see 
below however, is that some of the terms that are irrel- 
evant and neglected at the Luttinger liquid fixed point 
(which describes the paramagnetic phase only) are cru- 
cial for the stability of the ferromagnetic phase, as well 
as the ferromagnetic critical point; they must be retained 
for a proper description of the ferromagnetic phase as well 
as the transition. This should not be surprising, as oper- 
ators irrelevant at one fixed may well be relevant at other 
fixed points [12]. 

We now consider approaching the second order phase 
boundary from the paramagnetic side, by changing pa- 
rameters in the Hamiltonian. As we approach the criti- 
cal point, the spin susceptibility \ diverges; thus vn s cx 
l/x — * 0! As we move further into the ferromagnetic 
phase, one expects vn s to become negative; when this 
happens it becomes energetically favorable to have a non- 
zero expectation value of d x <j) s {x), which is the sponta- 
neous magnetization. Physically this occurs because the 
gain of exchange energy from the magnetization over- 
comes the loss of kinetic energy, as is standard in ferro- 
magnetism in itinerant electron systems. Clearly with a 
negative coefficient for (d x (j> s (x)) 2 , H s is not stable, and 
higher order terms in gradients of <j> s (x) and powers of 
d x (j) s (x) must be retained to maintain stability: 

H' = J dx[a{dlUx)) 2 + b(d x M*)) 4 + ■"], (?) 

where in the 

weak coupling limit a w (1/8tt 2 ) T,jf( V n,n+j ~ V n,n+ 3 ) 

and b w 2^^p^|fc=fc F ( e (k) is the single electron dis- 
persion). It is clear that the coefficient a is positive 
for generic repulsive interactions when there is Ising 
anisotropy (V A > V ++ = V~~). On the other hand 
the sign of coefficient b depends on the details of single 
electron dispersion; we assume b to be positive here [13], 
so that the transition is second order which corresponds 
to the situation in Ref. [5] . Other possible terms that are 
allowed by symmetry can be shown to be irrelevant at 
the ferromagnetic critical point to be studied below [14]. 
Thus combining H s with H', switching from Hamiltonian 
to action which is more convenient for RG analysis, and 
after rescaling the spin field (cf> — \/2a<f) s ), we arrive at 
the following effective action: 
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+ ^(d x cf>) 2 +u(d x cl > ) i 



(8) 



where r is imaginary time, r — v Ns /(ATta), and u = 
6/(4o 2 ). At mean field level, the ferromagnetic transi- 
tion occurs at r — 0. We note that the last three terms 
in Eq. (8) (or the static parts) take the usual form of 
Landau theory, upon identifying d x (f> as the local mag- 



netization m(x): 



i 4 . The first term 



in Eq. (8) controls the fluctuations along the imaginary 
time direction, and is responsible for the quantum na- 
ture of the theory. The effective action (8) is the basis 
of our study of the critical behavior of the ferromagnetic 
transition, which we now turn to. 

Renormalization Group and Critical Behavior 
at T = — We perform a renormalization group (RG) 
analysis of the action (8), and as is standard in such 
analysis, we treat the spatial dimension d as a contin- 
uous variable, even though the action (8) is derived in 
ID. Wc perform the following space-time transforma- 
tion with scale factor s > 1: x' = x/s,t' = t/s z , and 
(f)'(x') — (j)(x)s A , that leaves the first two terms in (8) 
invariant. This leads to z — 2 and A = d/2 — 1, which in 
turn lead to the (tree-level) scaling relation for r and u: 
r' = rs 2 and u' = us 2 ~ d . We thus find u is relevant be- 
low the upper critical dimension d c — 2, and the critical 
property of the transition is controlled by an interacting 
fixed point for d < 2. Obviously this fixed point is very 
different from the Luttinger liquid fixed point (which is 
non-interacting and has dynamical exponent z = 1), even 
though we are in ID. In the following we study the crit- 
ical property using momentum shell RG combined with 
e = 2 — d expansion, and hope it gives a reasonable de- 
scription at d = 1. 

We assume there is a momentum cutoff A, while no 
such cutoff exists in frequency. Integrating out modes 
with A/s < k < A and arbitrary frequency at one loop 
level, we obtain the following flow equations for r and u 
for small e: 
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dlogs 
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(9) 
(10) 



From these we obtain the fixed point (to order 0(e)): 
u* = (27r/9)e and r* = -(e/3)A 2 . This interacting fixed 
point controls the critical property of the transition, at 
which the RG dimension of the (relevant) tuning param- 
eter r is y r = 2(1 — e/6). From y r as well as the fact that 
the field 4> receives no anomalous dimension (77 = 0) to or- 
der O(e), we determine the critical exponents: correlation 
length exponent v = l/y r f=s (l/2)(l + e/6); susceptibility 
exponent 7 = (2 — n)v s=s 1 + e/6; magnetization expo- 
nent (3 = vd/2 w 1/2 — e/6; and field exponent 5 w 3 + e. 



To order 0(e 2 ) the field <j> receives a non-zero anoma- 
lous dimension 77, which also leads to a correction to the 
dynamical exponent z = 2 — 77/ 2 [15]. 

RG and Critical Behavior at Finite T — The 
phase transition occurs at T = only. However the 
quantum critical point has very significant influence on 
the thermodynamic properties at finite T, if the system 
is sufficiently close to the critical point. More specifi- 
cally, as we will see below, finite temperature introduces 
a thermal length scale £t ~ T~ 1/>z , and depending on its 
interplay with the correlation length £ of the system at 
T = 0, one can divide the temperature-coupling space 
into three regions, which are separated from each other 
by two crossover lines of the form T <~ \r — r*\ zv [16]: 

(i) Quantum Disordered: r > r* (where r* is the crit- 
ical coupling) and £t > £, in which the system behaves 
as an ordinary Luttinger liquid; for example the specific 
heat C ~ T and susceptibility \ ~ constant. 

(ii) Rcnormalized Classical: r < r* and £t > £, the 
system behaves like a ferromagnetic Luttinger liquid [8] 
whose magnetic order is suppressed by thermal fluctu- 
ations; here depending on the symmetry there can be 
two types of behavior: (a) if the system possesses full 
Heisenberg symmetry, the gapless transverse spin fluctu- 
ation with spectrum u> ~ k 2 gives rise to specific heat 
C - VT, while the susceptibility x ~ l/T 2 [17]; (b) if 
the system possesses Ising symmetry only, then the trans- 
verse spin fluctuation is gapped, the longitudinal fluctu- 
ation with linear spectrum (the usual Luttinger liquid 
behavior) gives rise to specific heat C <~ T, while the 
susceptibility \ ~ exp(J/T), which is the usual behavior 
of an Ising ferromagnet (J is an energy scale of order the 
domain wall energy of the Ising ferromagnet). 

(iii) Quantum Critical: £p < £, in which the thermo- 
dynamic property is controlled by the quantum critical 
fixed point studied above. In the following we focus on 
the quantum critical region. To study finite T proper- 
ties, we need to generalize the RG flow equation in the 
presence of finite temperature. The temperature T deter- 
mines the range of imaginary time: < t < l/T. Thus 
under scaling T scales as: T" = Ts z . The flow equation 
of r is modified to be 



dr 



2r^ A 2 - -r) coth — , 

ir y 2 ' 2T' 



dlogs 7r 2 ' 2T' ^ 

which leads to a flow away from the fixed point value r = 
r* due to finite T (flow of u is negligible if initially u = 
u*). Integrating Eq. (11) till r — r* <~ A 2 , at which scale 
the system is far away from criticality and the correlation 
length £(s) ~ 1/A, we find that the temperature T sets 
a correlation length 



as anticipated. This leads to the temperature depen- 
dence of susceptibility x ~ £ 2 ~ f / (eT 1 ) in the quantum 
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critical region, to order 0(e). To determine the behavior 
of specific heat, we need the singular contribution to the 
free energy from critical fluctuations, which obeys the 
hyper-scaling law for theories below their upper critical 
dimensions: 

p „T 1+d l*&(\r-r*\* v /T), (13) 

where <$>(x) is a universal scaling function. From Eq. 
(13) we immediately obtain C <~ T 1-e / 2 in the quantum 
critical region. 

Summary and Discussion — In this work we de- 
veloped a bosonic field theory that describes ferromag- 
netic transition in ID itinerant electron systems, based 
on Abelian bosonization. This approach is quite different 
from that of the Hertz-Millis theory, because in principle 
the bosonization procedure keeps all the degrees of free- 
dom of the original fermionic system, and allows us to 
arrive at a bosonic theory without integrating out gap- 
less fermionic degrees of freedom, as was done in the 
Hertz-Millis theory. We thus believe the theory devel- 
oped here is free of the possible singularities associated 
with integrating out gapless fermions. 

The bosonic theory developed here, when generalized 
to arbitrary dimensions, was found to have upper criti- 
cal dimension d c = 2. Thus for d = 1, which is where 
the theory applies, the system is below its upper critical 
dimension, and the universal critical behavior of the tran- 
sition is controlled by an interacting fixed point, which 
we have studied in some detail using e expansion. This is 
again quite different from the Hertz-Millis theory, which 
is above its upper critical dimension d c = 1 for d = 2 
and d = 3 where it applies; there the critical behavior 
is controlled by a Gaussian fixed point, and some of the 
critical properties are non-universal due to the presence 
of dangerously irrelevant operators. 

The Abelian bosonization procedure can be applied to 
systems either with full Heisenberg symmetry or Ising 
symmetry only. However since it does not exhibit the 
Heisenberg symmetry explicitly, this symmetry is most 
likely lost due to the approximate nature of the deriva- 
tion and treatment of the bosonic theory. Thus the crit- 
ical behavior discussed above probably applies to sys- 
tems with Ising symmetry only. In solid state systems 
the Heisenberg symmetry of electron spins are often re- 
duced to Ising due to the ubiquitous spin-orbit coupling; 
in systems with pseudospin transitions the Heisenberg 
symmetry is absent in the first place. Thus the results 
presented here are highly relevant. Nevertheless it would 
be highly desirable to maintain the Heisenberg symmetry 
when present, using for example non- Abelian bosoniza- 
tion, and study if and how the extra symmetry affects 
critical behavior of the transition. We leave this for fu- 
ture investigation. 

In an earlier work, Sachdev and Scnthil [18] studied 
ferromagnetic transitions in lattice rotor models, and ar- 
rived at an effective action very similar to Eq. (8); the 



only difference being the field is complex in their work, 
which is crucial for the Heisenberg (0(3)) symmetry of 
their model. If the symmetry were reduced to Z(2) xU(l), 
the corresponding theory would involve a real field, then 
the theory becomes identical to Eq. (8) [15]. It was con- 
jectured [15,18] that this action properly describes ferro- 
magnetic transition in ID itinerant electron systems. It 
is remarkable that two very different approaches lead to 
the same effective theory for the transition. 

Incidentally, a 2D version of the action (8) was used 
to study transitions between different valence bond solid 
states in 2D recently [19] . The physics of these transitions 
are very different from the one discussed here however. 

This work was initiated while the author was visiting 
the Max Planck Institute for Physics of Complex Systems 
in Dresden, during her Workshop on Quantum Phase 
Transitions. He benefited from stimulating discussions 
with Peter Kopietz. He is particularly grateful to Subir 
Sachdev and T. Senthil for explaining to him the connec- 
tion between Ref. [18] and the present work, and bringing 
Rcf. [19] to his attention. This work was supported by 
NSF grant No. DMR-0225698. 
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